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ABSTRACT 

A two-dimensional dynamic model of a towed buoy designed to remain 
submerged just below the water surface is developed. The study includes 
the model development, the analysis of certain model parameters , the 
development of a tension force, and the range of stability of the model to 
initial condition perturbations , 
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I. INTRODUCTION 



The purpose of this investigation is to attempt to mathematically 
model and computer-simulate the two-dimensional motion of a submarine- 
towed buoy. This has a very specific military application in submarine 
communications systems; the basic system configuration as a submarine 
communication link is illustrated in Fig. 1. It is desired that the buoy 
stay underwater at all times, and remain within ±1 ft. of four feet below 
the water surface; this gives the buoy an operating range from three to 
five feet below the surface. An additional system requirement is that the 
buoy system stay within its allowable operating zone while following 
waves up to sea-state four. 

Many investigators, such as Lacey, have dealt with ship-towed or 
airplane-towed devices, where the buoy is the "bottom-following" type. 
Pode and other early investigators studied the steady-state configuration 
of cables in a steady uniform stream; and more recently Schram developed, 
by the method of characteristics, a two-dimensional dynamic cable-buoy 
model. However, he treated the buoy as a vector force and did not go 
into the actual buoy dynamics . 

Steady-state studies do not show the response of the buoy and cable 
to perturbations such as sea state. It is this response that will prove 
most useful in the study of buoy and cable parameters. This study has 
been primarily concerned with developing the dynamic model of the buoy; 



and the tension force generated by the cable has been loosely approximated. 
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SUBMARINE COMMUNICATIONS BUOY 



The model is of a buoy presently being developed for the U. S. Navy; 
thus , the dimensions and other data were obtained from an existing 
prototype . 

This study develops the two-dimensional model around two basic 
equations; force equals mass times acceleration, and the moment equals 
the mass moment of inertia times the angular acceleration. The model 
block diagram was built around these two equations, several basic 
hydrodynamic relationships , and the wind-tunnel and dimensional data 
supplied on the buoy. Once the block diagram was formulated, it was 
then ready to be computerized. The model was simulated using the IBM 
360-67 at the Naval Postgraduate School, and DSL (Digital Simulation 
Language); DSL is a very simple, but versatile simulation language 
developed by IBM. Once the model was computerized, the lengthy and 
time-consuming task of optimizing certain parameters and finding a real- 
istic formulation of the tension force followed. The determination of the 
tension force and its coupling with the buoy has been one of the most 
difficult problems in the formulation of the buoy-cable model. This has 
been a problem in all buoy-cable models where the dynamic response is 
of interest; this is mainly due to the fact that it is not exactly known how 
the forces generated by the buoy and cable react on and with each other. 
The development of the model and the results of the simulation are 
presented . 
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II. EQUATIONS OF MOTION AND PROBLEM SETUP 



The purpose of the buoy is to remain four feet below the water 
surface so that a protruding antenna may maintain two-way communica- 
tions; when the buoy is operating in high seas it must follow the sea 
surface. The body, as shown in Fig. 2. , is basically a fat wing that 
produces dynamic lift according to its angle of attack. This angle of 
attack may be changed by torque produced by the fin. Besides the 
dynamic lift, produced by the body angle of attack, there is the static 
lift of positive buoyancy. The buoy model is based on the two basic 
equations ; 

F=ma 

M=I V 
y 

The first equation is used to describe movement in the horizontal and 
vertical planes, the X and 2 planes respectively. The second is used 
to describe the rotation of the body about the Y axis; these axis are 
shown in Fig. 2. The three hydrodynamic forces of lift, drag, and moment 
are defined as follows; 

L= 

D=iCj^pV"S 

M=Kj^pV"SC 

Along with the tension force produced by the cable and the buoyancy of 
the body, these three equations describe all the forces acting on the 
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THE BUOY 



body. Before going into the complete buoy block diagram, the body 
coefficients will be put into a more convenient form, and the fin servo 
loop, which controls the buoy attitude, will be examined. 

A. BODY COEFFICIENTS 

In both aerodynamics and hydrodynamics , there are various dynamic 
coefficients; these are determined by the shape of the body involved. 

There are tables of coefficients for certain common wing shapes. However, 
when an odd body configuration is being studied, it may not be possible 
to apply the standard shapes in order to derive overall coefficients for 
the body. In this case the coefficients of the body must be determined 
by wind-tunnel tests. 

This was the case with the system being studied; the body is of 
non-standard shape, but the fin is the standard NACA 0012 shape. The 
three coefficients, C , C , and 0 , were measured at various body 

i-/ IVl 

and fin angles. The coefficient of lift and the coefficient of moment 
include the effect of the positive buoyancy force of 490 pounds. The 
effect of the antenna drag is not included in any of the coefficients . 

The values of these coefficients are shown in Tables I-III as a function 
of body angle of attack and fin angle of attack. It was decided that for 
computer implementation, these coefficients would be easier to handle 
if they were in equation form. This was done by applying a least-squares 
fit twice to each data set; the resulting equations may be found in Sub- 



routine Tink of the digital program of the model. The resultant curves and 
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data points for constant body angle of attack and varying fin angle of 



attack are shown in Figs. 3-5. 

B. FIN SERVO LOOP 

The diagram of the fin and the servo-control system is shown in 
Fig. 6. The purpose of the fin is to orient the body so that it will stay 
four feet below the water surface. The fine produces a moment on the 
buoy which causes the body angle of attack to increase or decrease. 
This increase or decrease in body angle of attack causes the body lift 
to increase or decrease and the body rises or descends as is necessary. 

The fin actuator loop is driven by a Size 23, Kearfott motor. The 
loop is designed to have a rise time of 0.04 sec.; this actuator loop is 
shown in Fig. 7. The servo loop used in the buoy model is considered 
to be a K/s^ plant with tachometer and unity feedback; this is shown in 
Fig. 8. The system transfer function: 



G 




1 + GH 



1 + K^/s=(KjS + 1) 



s= + 



and the characteristic equation; 



s^ + K„K s + K = 0 . 

A t A 



Assuming K =1000 

/A 



w^ = K = 1000 
n A 
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the fin servo control system Angle Sensor 

Figure 6 





THE FIN SERVO LOOP 
Figure 7 




THE K/S^ SERVO SYSTEM 
Figure 8 
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24w = K„K 

n At 



K 



t 



24 w 

n 

1000 



Taking a value of 4=0.7: 

= 2(0.7)(31.62)/1000 

= .0443 . 



With K- and K set at the above values , the rise time of the servo 
A t 

loop is about 0.1 sec. This is not as fast as desired, but by changing 
the two gains the desired rise time may be obtained. Two different 
responses of the fin servo loop are shown in Fig. 9; both have an input 
step of 0 . 1 74533 rad . (10 degrees). The first set of gains is used in the 
model simulation. This is because it is doubtful that the inputs into the 
servo loop will require that the fin move that fast; if so, then the buoy 
has probably reached the limits of its surface-following capability. 

The inputs into the fin servo loop are shown in Fig. 6. The fin 
servo loop that was simulated above is composed of the summer, servo 
amplifier, motor and gears, linkage to the fin, and the unity and tacho- 
meter feedback from the fin to the summer. The other three inputs will 
come from gages and electronic sensors in the body. These inputs are 
the depth error, the depth rate, and the body pitch rate. In the model 
simulation, there arc three constants (ADE, ADR, and APR) which multiply 
these three inputs. The values of these three constants are discussed 
later in the study. 
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Figure 9 



27 



III. THE BUOY MODEL 



In this development of the buoy model, the effect of water currents, 
the forces generated by the heaving water surface, and the free-surface 
effect are not taken into account. Currents can be considered neglible 
compared to the heaving forces the buoy might experience. The free- 
surface effect is experienced by the buoy because it is a body moving at 
a depth comparable to its length; a more detailed and complicated model 
would have to include this force. The force generated by the heaving 
water surface would not affect the buoy as much as might be expected. 

This is because the body, as lonq as its vertical freedom is not restricted, 
will act like a particle in the wafer and will move up and down with the 
surface. Thus, the heaving forces will actually help maintain the buoy 
near the water surface, which is what is desired. If the buoy were being 
towed directly above its towpoint, it would have little or no vertical 
freedom; however, if the buoy is well behind the towpoint, then the cable 
will provide little restriction to the buoy's vertical freedom. Tests on the 
actual buoy have shown that with 475 ft. of cable and the two point at a 
depth of 240 ft. , the buoy can heave up to 20 ft. with the waves without 
the cable decreasing the buoy's vertical freedom of motion by more than 
10 percent.^ 



^ Soulant, H., Sea Surface Followaqc by Towed Buoy (Preliminary 
Study) , p. 11, March 19 69. 
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The buoy model will now be developed; the buoy diagram is again 
produced with some specific dimensions (Fig. 10). Besides the diagrammed 
dimensions, the buoy has the following additional characteristics: 

S = 14.25 ft^ 

C = 6 . 0 f t 

b = 2.375 ft 

mass = 70 slugs 

mass moment of inertia = 975 slug -ft® 
volume = 17.5 ft® 

The block diagram of the dynamic model is shown in Fig. 11. As 
was previously stated, the buoy model is based on the summation of 
forces in the X and Z planes, and about the Y axis . The summed forces 
in the X and Z directions are divided by the mass of the buoy, and the 
moments about the Y axis are summed and divided by the mass moment 
of inertia. This yields the acceleration in the horizontal and vertical 
directions , and the angular acceleration about the Y axis . The accelera- 
tions are each integrated to give the velocities, and these, when inte- 
grated, give the horizontal and vertical position and the angular 
orientation. 

In this first model, the tension force is generated as a function of 
horizontal velocity and body angle of attack. The tension function was 
formed in the same manner as the coefficient equations , and was based 
on predicted values of the tension in the cable. This tension function may 
also be found in Subroutine Tink of the computer program. 
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THE MODEL BLOCK DIAGRAM 
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THE MODEL BLOCK DIAGRAM 





The velocities and positions are used to generate the forces and 



moments that are summed in the basic equations. The forces and moments 
on a body are relative to the velocity vector of the medium , and not to the 
local horizontal. Thus, the sum of the squares of the horizontal and 
vertical velocities yields the relative velocity squared, V2 . The arc 
tangent of the vertical velocity over the horizontal velocity is taken to 
give the, angle of the velocity vector, VAN, with respect to the true 
horizontal. The angular orientation of the body (angle of attack), A1 , is 
summed with VAN to produce the angle of attack of the body with respect 
to the velocity vector, AR; AR is in radians. The fin angle of attack, FA, 
is also summed with VAN to yield the fin angle of attack relative to the 
velocity vector, BR. These two angles of attack, AR and BR, are 
multiplied by a conversion factor of 57.29578 deg/rad to yield the angles 
of attack , A and B , in degrees . 

The two angles, A and B, are the arguments which determine the 
coefficient values of body lift, drag, and moment. The coefficients of 
lift and drag are both multiplied times V2 and a constant to produce the 
body lift, BL, and the body drag, BD. The constant is ^pS, and is 
equal to 14.17875. 

The moment coefficient is multiplied times V2 and a constant, •^pSC; 
this constant is equal to 84.0725. This body moment, YCM, is summed 
with the moments produced by the antenna, to yield the total moment, YM, 
acting on the body about the towpoint. 



.3 



There is no definite information on the antenna, but it is assumed 



that it will only contribute a drag force; this drag force will contribute to 
the total drag and moment on the body . The antenna coefficient of drag is 
assumed to remain constant at a value of 0.50. Thus, the antenna drag 
is dependent on a constant times the frontal surface area times the 
relative velocity squared. Assuming a frontal width of 2 inches, the 
drag force of the antenna will then vary with the length of antenna under 
water, and V2. The length of antenna will be approximated as the depth 
of the body below the surface. The constant in this case is ^(l/6)pC , 

J— 

which equals 0.04146. The antenna drag is multiplied by the sine and 
cosine of VAN, SV, and CV, to give two moment components. These times 
their respective moment arms yield YADS and YADC . 

The body lift and drag, and antenna drag are multiplied by SV and CV, 
and the resulting six values are summed into the horizontal and vertical 
forces. In this first realization of the model, these horizontal and 
vertical forces are used before being summed with the horizontal and 
vertical components of the tension force, to determine the angle FE with 
which the tension force acts on the body. This means that the tension 
force acts exactly opposite to the vector of body forces; this simplifying 
assumption is used in the initial stages of the model development. 

The depth error and depth rate inputs into the fin servo loop are taken 
from the vertical position and velocity, Zand Zl , respectively. These 
arc multiplied by a constant, 0.443, which converts ft and ft/sec into 
psi/ft and psi sec/ft. These arc multiplied by gain constants ADE and 
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ADR before being summed with the body pitch rate, Y1 , times its gain 



constant, APR. These three inputs are fed into the fin servo loop. The 
orientation of the fin causes a torque on the body which either increases 
or decreases the body angle of attack; this enables the body to rise or 
descend. 

This is the initial dynamic model; it changes as more insight and 
knowledge are gained on the model. The first study with the model 



deals with the three coefficients ADE, ADR, and APR. 



IV. STUDY OF THE BUOY MODEL 



Once the model was digitally programmed and the initial programming 
mistakes ironed out, the first step was to study possible values of the fin 
servo input coefficients. The first computer program realization is shown 
in the program titled "INITIAL BUOY MODEL" . The buoy is initiated at a 
depth, ZI, of 4.0 ft, a speed, XII of 20.2536 ft/sec (12 knots), body 
angle of attack. All, of 0.10472 rad (6 degrees), and a fin angle of attack, 
FAI, of 0.02 rad (1.15 degrees). The tension force and the coefficients 
of lift, drag, and moment are called out of Subroutine Tink, which is in 
double-precision format. 

A. SERVO INPUT COEFFICIENTS STUDY 

The first attempt in this study was to see if the buoy would stabilize 
with only one input into the fin servo loop; either the depth error, the 
depth rate, or the body pitch rate. This proved fruitless, so different 
combinations of ADE and ADR were tried, and it was found that these two, 
with APR set equal to zero, would stabilize the system. A full range of 
values of ADE and ADR were tried to see which combinations would bring 
the buoy to a steady-state position. A mapping of these values is shown 
in Table IV; the value given for each combination is the depth at which 
the buoy steadied out. The responses for several of these combinations 
are shown in Figs. 12-16. It appears from the depth values that the 
tension force is too strong, since in every case the buoy steadies out 
deeper than the desired depth of 4.0 ft. This became obvious when the 
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MAPPING OF STABLE COMBINATIONS OF ADE AND ADR 
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RESPONSE OF ADE=0.50 AND ADR=0. 
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RESPONSE OF ADE=1.6 AND ADR=0,09 
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RESPONSE OF ADE=1.2 AND ADR=0 
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RESPONSE OF ADE=1.2 AND ADR=0.07 
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RESPONSE OF ADE=2.5 AND ADR=0 



steady-state body angle of attack and fin angle of attack were checked; 



both angles were large and showed that the model was working properly, 
but that the buoy was fighting against an unrealistic tension force 
formulation. 

It was decided from the mapping of the two coefficients, ADE and 
ADR, that three sets would be used in further studies of the buoy model. 
The first pair was chosen on the basis that it was on the intersection of 
the values of ADE and ADR which had the widest range of stable values , 
and that this set was towards the center of the mapping. This was 
ADE-1 .2 and ADR=0.30. It was then decided that the other two pairs 
would vary in one coefficient away from this point. These other two 
pairs are ADE=1.2 and ADR=0.07, and ADE=2.5 and ADR=0.30. The 
responses with these three pairs of coefficients are shown in Figs. 

14-16. 

The coefficient APR remained to be studied. When this coefficient 
was given a value larger than zero there was no appreciable change in 
the response of the buoy. To check the range of possible values of APR, 
the tension force, T, was changed so that it was equal to the square 
root of the sum of the squares of Fi and F2 . Thus, the buoy would stay 
at four feet below the surface and the effect of APR could be seen in the 
response of the buoy's angle of attack. Both negative and positive values 
of APR were tried, and negative values immediately drove the buoy model 
unstable. Positive values of APR, up to 7.0, were tried with the three 
pairs of depth error and depth rate coefficients . At this upper value of 



4 3 



7.0, the body started to oscillate away from stability about the Y axis; 
it was apparent that this high value of gain was approaching the stability 
limit of the buoy model. Based on this, a median value of 3.0 was 
chosen for APR. The three final sets of coefficients are as follows: 



ADE 


ADR 


APR 


1.2 


0.07 


3.0 


1.2 


0.30 


3.0 


2.5 


0.30 


3.0 



B. TENSION FORCE FORMULATION 

The next step in this model study was to formulate a new, and hope- 
fully more realistic, tension force, rather than using the predicted 
tension data supplied with the buoy. It was decided that this tension 
function would be a constant plus or minus the distance away from the 
equilibrium point times what will be called a spring constant, SC. The 
basis for this formulation is shown in Fig. 17, First it is assumed that 
the last portion of the cable connected to the buoy, LO, is straight, and 
that of the total length of the cable this section is the most dynamic. By 
computing the increase or decrease in ZO and XO, a new cable length, 

LN, may be calculated. If LN is larger than LO then the tension increases; 
if it is smaller then the tension decreases. The difference in length 
(LN-LO) is L and this is multiplied times the spring constant, SC, which 
is studied for several different values later on. The new tension force 
now takes the form: 

T = 1940 + L-SC . 
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THE NEW TENSION EORCE EORMULATION 
Figure 17 
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The value 1940 is the steady-state tension value that the model 
settled to in the earlier parameter studies; the angle FSS is the steady- 
state value of FE, This new tension force was implemented into the buoy 
model. The buoy will now be initiated at a depth of zero feet; this is 
possible because no heaving force or surface effect is taken into account 
in the model; other minor changes are made to take this new initial 
condition into account. The new computerized model program is titled 
"STUDY OF T = 1940 + L*SC". The body angle is now being initiated at 
4.8 degrees and the initial fin angle is 0.0 degrees; these two values 
are also the steady-state values for these angles found earlier. The next 
step in the formulation of this tension force was to find an appropriate 
value for the spring constant. The program was run with the three sets 
of coefficients and the new initial conditions , while SC was varied from 
2000. to 0.0 in steps of 250. The two sets of coefficients with ADE=1.2 
were stable for all values down to SC=250, while the set with ADE=2.5 
was stable down to SC=1750. This third set (ADE=2 . 5 , ADR=0.30, and 
APR=3.0) was dropped from further investigation since this was beUeved 
to be an unrealistic value for the spring constant. 

C . STUDY OF INITAL CONDITION PERTURBATIONS 

This final form of the model with the new tension force formulation 
was next examined to see how stable the model was under various initial 
condition perturbations. It was decided to perturb the initial vertical 
velocity, the body angle of attack, and the fin angle of attack. Three 
diflcrcnt spring constants were used with the two lemaininy sets of 
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coefficients; these were SC=500, SC=1000, and SC=1500. It was also 



decided that two different values for the approximation of the dynamic 
length of the cable would be used; these were ZO=100 and ZO=50. (ZO 
simply because it was used to initiate the model rather than LO) . The 
final results were tabulated, and are shown in Table V; the initial angles 
of attack (All and FAI) are in degrees , and the initial vertical velocity 
(Zll) is in ft/sec . 

The responses of the buoy model to some of these initial condition 
perturbations are shown in Appendix A. The twenty-four plots include 
all the maximum perturbations listed for SC=500 and SC=1500. These 
plots do not really tell much about the buoy model. They do show that 
the buoy will settle to steady state fairly rapidly, and that the buoy will 
not deviate very far from its desired depth. Plots with the same pair of 
coefficients (ADE and ADR), with the same dynamic length of cable, and 
with the same type of inital condition perturbation, but with different 
values of spring constant, seem to be very similar. The responses with 
the larger spring constant (SC=1500) seem to settle more quickly, even 
though they are oscillating through a wider range of depths . All the 
responses to the fin angle perturbations are remarkably similar. This 
might be because the fin is quickly driven back to the desired stability 
position; and since all the fin angle perturbations are the same pair of 
values, this similarity in responses is not surprising. The perturbing of 
the initial conditions for the new tension force realization was the point 
at which this study ended. 



47 



ZO=50. 



Z0=100 



SC=500 


ADE=1 .2 
ADR=0.07 
APR=3 . 


ADE=1 .2 
ADR=0 . 3 
APR=3 . 


ADE=1 .2 
ADR=0.07 
APR=3 . 


ADE=1 .2 
ADR=3 . 0 
APR=3 . 


All 


-2. 1-+0.9 


-1 .8-+0.9 


-1 .7-+1. 


-1 .5-+0.9 


All 


+2 . 5-*+8 . 0 


+3 . 0 — +7 . 5 


+2 . 0— +7 . 


+2 . 5 — +6 . 5 


FAI 


-17.5-H5.0 


-17 .5-+15.0 


-17 .5-+15.0 


-17.5-+15.0 


SC=1000 


Zll 


-4.8-+2.S 


-2 .9-+2 .6 


-4.3-+2.5 


“2 , 8“^ +2 . 2 


All 


-0.5-+9 .5 


0 „ 0— +9 . 0 


—0 . 5 — +8 . 0 


0 . 0“^ +8 . 0 


FAI 


-17.5-15.0 


-17 .5-+15. 0 


-17.5-+15.0 


-17.5-+15.0 


SC=1500 


Zll 


-4 . 8-* +3 . 1 


-2 .7-+3 .3 


-4.8-+3.0 


-2 .7-+2 .9 


All 


-2 .0- + 10.0 


-1 .0-+9.5 


-1.5- +9.0 


-1.0-+8.5 


FAI 


-17.5-+15 .0 


-17.5-+15 .0 


-17.5-+15.0 


-17.5- + 15.0 



STABILITY RANGE OF INITIAL CONDITIONS 
Table V 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The study presented here is only a beginning to this problem. There 
are many facets and additions to this problem that remain to be studied. 
First of all, the dynamic characteristics and wave-following capabilities 
of this model were never investigated. Other additions could be made to 
the model to take into account the heaving forces generated by the waves 
and the free-surface effect which will definitely affect the motion and 
capabilities of the buoy. There remains much to be studied in the area 
of the dynamic cable and the coupling of the tension force to the buoy 
model. Other tension force realizations might be tried with several spring 
constants in a single cable length; or the cable might be approximated by 
a series of straight lengths, each with its own spring constant. The 
possibilities of adding a damping coefficient into the tension force real- 
ization need to be investigated; the possibilities are almost endles . 

Other methods of simulating the buoy need to be tried, especially that 
of using a hybrid computer (an analog computer integrated with a digital 
computer); hybrid computers are more suited to this type of real-time 
simulation work. A dynamic model that allowed for all the possible 
dynamic forces and all the possible variations in buoys and cables, 
would be invaluable in the design and engineering of towed buoys . 
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APPENDIX A 



RESPONSES OF THE BUOY TO INITIAL CONDITION PERTURBATIONS 
The following 24 plots are the responses of the buoy to the maximum 
initial condition perturbations listed in Table V for SC=500 and SC=1500. 
The plots are of depth versus time. The x scale for all the plots is 
5 sec/inch; and the y scales are in ft/inch. 
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APPENDIX B 



SYMBOLS AND ABBREVIATIONS USED IN THE COMPUTER PROGRAMS 
Z vertical position in ft. (initial condition ZI) 

Zl vertical velocity in ft/sec. (initial condition Zll) 

Z2 vertical acceleration in ft/sec^ 

FV summation of vertical forces 

FI summation of the vertical forces less the vertical tension 

com ponent 

XI horizontal velocity in ft/sec. (initial condition Xll) 

X2 horizontal acceleration in ft/sec^ 

FH summation of the horizontal forces 

F2 summation of the horizontal forces less the horizontal 

tension force 

FF the absolute value of FI 

FG the absolute value of F2 

FE the angle with which the tension force acts on the body 

in radians 

SFE the sine of FE 

CFE the cosine of FE 

TSF the vertical component of the tension force 

TCF the horizontal component of the tension force 
C horizontal velocity in knots 

VAN the angle of the velocity vector relative to the horizontal 
in radians 

V2 the magnitude of the velocity vector squared 
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A1 

Y1 

Y2 

YM 

AR 

A 

FA 

FB 

FC 

FD 

BR 

B 

CL 

CD 

CM 

T 

CV 

SV 

ZAB 

AD 



the body angle of attack relative to the horizontal in radians 
(initial condition All) 

the body pitch rate about the Y axis in radians/sec 
(initial condition Yll) 

angular acceleration of the body about the Y axis in 
radians/second^ 

summation of the moments on the body about the Y axis 

the body angle of attack relative to the velocity vector 
in radians 

AR in degrees 

the fin angle of attack relative to the horizontal in 
radians (initial condition FAI) 

used in the fin servo loop 

used in the fin servo loop 

summation of the inputs and the feedback in the fin 
servo loop 

the fin angle of attack relative to the velocity vector in 
radians 

BR in degrees 

coefficient of lift of the body 

coefficient of drag of the body 

coefficient of moment of the body 

the tension force 

the cosine of VAN 

the sine of VAN 

the absolute value of z 

the antenna drag 
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BL the body lift 

BD the body drag 

YCM the moment on the body about the Y axis 

YADS the moment on the body due to the vertical component of 
the antenna drag 

YADC the moment on the body due to the horizontal component 
of the antenna drag 

ZD in the first program this is used to measure the relative 

depth error from 4,0 ft; in the second program it is used to 
compute the antenna drag 

DEE the depth error in psi/ft 

DE DEE times the coefficient ADE 

DRR the depth rate in psi* sec/ft 

DR DRR times the coefficient ADR 

YPR the body pitch rate (Yl) times the coefficient APR 
ADE the depth error coefficient 

ADR the depth rate coefficient 

APR the pitch rate coefficient 

THE FOLLOWING ARE USED SOLELY IN THE SECOND COMPUTER PROGRAM 
IN THE REALIZATION OF THE NEW TENSION FORCE 

FSS the steady-state value of FE 

ZO the vertical approximation to the dynamic length of the 

cable; see Fig , 1 7 

SC the spring -constant value of the cable in pounds/ft 

SFS the sine of FSS 

LO the approximation to the dynamic length of the cable; 

see Fig . 1 7 



TFS 

XO 

ZN 

XR 

X 

XN 

LN 

L 

FE 



the tangent of FSS 

the horizontal approximation to the dynamic length of the 
cable; see Fig . 17 

the new computed value of ZO 

the relative movement out of steady state in the X 
direction 

the change in position out of the steady state in the 
X direction 

the new computed value of XO 

the new computed length of the dynamic length of the 
cable 

the change in length of the cable 

as stated previously, but computed differently 
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CCMPUTATICN OF ANGLES OF ATTACK RELATIVE TO THE VELOCITY VECTOR 

AR=A1+VAN 

A=AR*57. 29578 

BR=FA+VAN 

B=BP»57. 29578 
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NOTE — BLANK CARDS AND CARDS MARKED ARE NOT INCLUDED IN ACTUAL PROGRAM 
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